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Abstract 

We calculate certain features of Bose-Einstein condensation in the ideal gas by using recur- 
rence relations for the partition function. The grand canonical ensemble gives inaccurate results 
for certain properties of the condensate that are accurately provided by the canonical ensemble. 
Calculations in the latter can be made tractable for finite systems by means of the recurrence 
relations. The ideal one-dimensional harmonic Bose gas provides a particularly simple and peda- 
gogically useful model for which detailed results are easily derived. An analysis of the Bose system 
via permutation cycles yields insight into the physical meaning of the recurrence relations. 
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I. INTRODUCTION 



The achievement of Bose-Einstein condensation (BEC) in alkali gaseail was a remarkable 
feat in atomic and low-temperature physics. The gases are most often trapped magnetically 
in potentials accurately approximated by harmonic oscillator wells. The result has been a 
deluge of theoretical papers on BEC in harmonic potentials, both for ideal and interacting 
gases.i Experimentally, the gases are sufficiently dilute and weakly interacting that the ideal 
gas is a good first approximation for their description, making the subject more accessible 
to students even in their first statistical-physics course. 

A basic problem with the standard presentation of BEC is that the grand canonical en- 
semble misrepresents several physical quatities when a condensate is present. For example, 
grand canonical ensemble calculations greatly overestimate the fluctuations of the conden- 
sate number. Although, because of their isolation, the most realistic description of the 
experimentally condensed gases is via the microcanonical ensemble, the canonical ensemble 
gives equally accurate results. In the latter the system of interest is in contact with a heat 
bath but the particle number is kept fixed, which is crucial. 

Calculations using the canonical ensemble are avoided in most elementary treatments 
of BEC because of mathematical complications. The grand canonical ensemble removes 
these complications by putting the system in contact with a particle bath. Unfortunately, 
when there is a condensate, the de Broglie wavelength can be larger than the system size, 
making the distinction between system and bath meaningless and leading to the fluctuation 
inaccuracy mentioned. Moreover, the trapped Bose systems of current interest consist of 
relatively few particles. There exist simple methods using recurrence relations that have been 
exploited often in recent work to treat finite systems in the canonical ensemble; however, 
these relations are not well known outside the research literature. It is these techniques that 
we want to discuss here. 

Because the experimental trapped systems are finite, they are not equivalent to systems in 
the thermodynamic limit. The trapped systems have a nonuniform distribution established 
by the external harmonic potential. The true thermodynamic limit in these systems would 
require that as the particle number N is increased, we also would decrease the frequency uo of 
the potential in such a way that the maximum density — proportional to u d N, where d is the 
spatial dimension of the system — remains constant Although this limit is not physically 
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realized and finite systems have no real phase transitions, the experimental transformation 
of the system into its lowest state is still rather sudden. Nevertheless, the existence of 
a mathematically sharp phase transition is not crucial to the description of real systems. 
What is important is the appearance of a "condensation," by which we mean the rapid 
accumulation of a substantial fraction of the iV particles into the ground state (without big 
fluctuations about this average) when the temperature falls below a certain finite value. We 
will show that even a finite one-dimensional ideal Bose gas in a harmonic potential has this 
property. 

The choice of the canonical ensemble and the use of the recurrence relations are particu- 
larly suitable for the study of finite systems. Thus a simple model with physical properties 
amenable to calculation is available to be exploited for pedagogical or other purposes. 

In Sec. |n] we will review the standard grand canonical ensemble treatment of the one- 
dimensional (ID) harmonic ideal Bose gas and identify a temperature below which there 
is a substantial accumulation of particles in the ground state. We will also identify the 
physically unrealistic fluctuations in the ground-state occupation that appear in the grand 
canonical ensemble description. In Sec. [TTI| we will show how the ID Bose gas can be treated 
by developing a recurrence relation for the partition function. More general recurrence 
relations for the average number of particles in a single-particle state and for the partition 
function are developed in Sec. [TV]. Armed with these tools, we compare canonical ensemble 
calculations with their grand canonical ensemble equivalents. For a finite system there are 
only small differences in the mean values; however, there are large differences in the root- 
mean-square fluctuations. In Sec. [V] we look at the condensation problem from a quite 
different point of view, namely that of permutation cycles a la Feynman.! Such a view gives 
a physical explanation to several mathematical formulations found in the previous sections 
and especially to the partition function recurrence relation. From this point of view we also 
see that the grand canonical ensemble misrepresents the condensate, while the canonical 
ensemble treats it accurately. We find the somewhat surprising result that the condensate is 
made up of equally probable permutation cycles of all lengths up to the condensate number. 
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II. GRAND CANONICAL TREATMENT 



We first consider a grand canonical ensemble treatment of BEC in a one-dimensional 
(ID) harmonic well. Although our approach is typical of most statistical physics textbooks, 
we know of only one such book that actually covers this particular example.! The harmonic 
potential leads to equally spaced single-particle energy levels given by 

e p = pA, (1) 

with p a nonnegative integer. The zero-point energy, omitted in Eq. (|]), can be restored 
to any physical quantity at the end of the calculation. The constant A is related to the 
harmonic angular frequency u by A = Uu. 

In the grand canonical ensemble! the average number of particles N is given by the 
relation 

p 

where (3 = 1/ksT and /x is the chemical potential. In Bose problems the denominator is 
often expanded in powers of e-^-^ to yield 

oo oo oo 

N = £ e 1 ^ £ e"^ A = £ e l ^Z l ((3l), (3) 
1=1 p =o 1=1 



where 



z m = t^t^a ( 4 ) 



is the one-body partition function at the effective inverse temperature (31. In Sec. V we will 
see that the sum over I in Eq. (|3|) represents a sum over permutation cycles. 

For the very weak potentials used to trap iV particles experimentally, the harmonic os- 
cillator states are very closely spaced (hu <C ksT). Thus to a good approximation we can 
replace the sum over p in Eq. (|]) by an integral. (Alternatively, we could directly replace the 
sum in Eq. (pj) by an integral to arrive at the same result.) Because Jdpe~^ pA = (/3ZA)~ X , 
we find 

1 oo i b T 

N ~ N > = -Ly e W± = i n (i _ g/JM). ( 5 ) 

(3 A jr[ I A 1 w 

The sum in Eq. (|5]) is one of the Bose integrals,! which in this case can be evaluated 

analytically. Of course, changing the sum to an integral is valid only if the summand is 

a smooth function of p; we then lose the contribution of the lowest state when it becomes 



occupied with order N particles. Hence N' in Eq. fl5|) is just the contribution of the excited 
states, and we obtain the total number N by including the ground-state population: 

N = n + N', (6) 

where 

n = (e-^ - l)" 1 . (7) 

We want to identify a Bose-Einstein "transition temperature," that is, one below which 
there will be a sizable fraction of the N particles in the ground state. From Eq. (|7]) this 
requirement implies that —/3fj, = ( , ~yN)~ 1 <C 1, where 7 is a number of order unity. If we 
invert Eq. we obtain 

1 - e-? N ' A = « 1 + Pfi, (8) 

so that /3N'A pa In^N. Because N' is of order N at the transition, we find that the 
condensate will be large for temperatures below T defined byi! 

We can show! that the density of the system is proportional to iVA, so that, in the ther- 
modynamic limit, we keep the numerator iVA constant while letting N — > 00, A — > 0. Then 
the characteristic temperature will go to zero as 1/lniV, which is small only for extremely 
large N. In actual ID experiments,0'0 where N is about 10 4 , the logarithm reduces this 
characteristic temperature only by a factor of order ten compared to iVA/fc^, an experimen- 
tally accessible value. Nevertheless, we say that there is a "quasi-condensation" rather than 
a real one. For a two-dimensional (2D) ideal gas we find that an actual phase transition 
occurs; however, in accord with the Hohenberg theorem,! this transition disappears if there 
are particle-particle interactions. On the other hand, some authors have claimed that the ID 
and 2D finite interacting systems at sufficiently low temperature have a true condensation, 
because the coherence length becomes larger than the finite condensate sizeJllS (Moreover, 
in 2D we expect a true phase transition of the Kosterlitz-Thouless varietyB at a temperature 
of order As/N.) 

It is straightforward to solve Eq. (|3|) numerically for the chemical potential fi and compute 
the exact condensate number Eq. fl7|) for some given average iV value. In Fig. 1 we show the 
occupation of the first two energy levels for iV = 500. We see that To provides a fairly good 
estimate of the quasi-condensation temperature. 



Many textbooks compute particle-number fluctuations in the grand canonical ensemble.! 
The grand partition function for any ideal Bose gas with states e p each occupied by n p 
particles i si 

OO -| 

z = Y,z N = n E = n T—^ry ( 10 ) 

N V n p =0 p L e 

where is the canonical partition function of an iV-particle system. The average square 
deviation of the occupation number n p is 

A^2 = (n p -n p y = (k B T) 2 ^^ = n p + n\. (11) 

(Jt p 

With no condensate we have n„ C 1 for all except a negligible number of excited states, so 



AN 2 = J2 P An p ~ Yl p n p — N, and we have a normal distribution with 



N 

However, with a condensate of order N, we have 



^ =0 (N-^). (12) 



-^ = 0(1). (13) 

that is, the fluctuations of the condensate are as large as the condensate itself — a manifestly 
unphysical result. This problem is not new, but has received a large amount of recent 
attention00@@0000ii@@ including the invention of a new "fourth" ensemble— the 
so-called "Maxwell demon" ensemble — to take care of it.0 

There are various explanations of what goes wrong with the grand canonical ensemble. 
Grossmann and Holthausll state that "[T]he relative mean square fluctuations of the ground 
state population, and thus the relative fluctuations of the total particle number, approach 
unity: as a result of particle exchange with the reservoir, the uncertainty of the number 
of particles becomes comparable with (N) itself. This fluctuation catastrophe is related to 
the divergency of the quantum coherence length At for T — > 0. When At vastly exceeds 
the length scale characterizing the system under consideration, a rigid distinction between 
'system' and 'reservoir' is no longer practical." The difficulty also can be stated in more 
mathematical terms. The grand canonical ensemble often is shown to be equivalent to the 
canonical ensemble by using the method of steepest descents,!! which evaluates the canon- 
ical partition function by an approximation to a complex integral of the grand canonical 
partition function. Wilkens and Weiss^l state that "The most common procedure is to eval- 
uate the contour integral in a stationary phase approximation. In leading order one recovers 



the grand-canonical formulation. However, below T c , fluctuations are badly represented in 
this approach. The reason is that, for large N, the saddle point is located within a dis- 
tance 0(1/N) from the branch point while the Gaussian approximation for the fluctuations 
assumes a much larger range of validity 0(1/ 'iV 1 / 2 )." 

We will not pursue the cause of this difficulty further, but will avoid it by using the 
canonical ensemble. 



III. ONE-DIMENSIONAL BOSE GAS BY THE CANONICAL ENSEMBLE 



We next examine the ID ideal Bose gas by using the canonical ensemble. There is a 



simple recurrence relation for the partition function of this model. As we will see in Sec. |TV| , 
there are more general recurrence relations by which other problems (for example, the 3D 
ideal Bose gas) can be treated. The connection between the two formulations turns out to 
be a special case of a famous theorem in number theory as we demonstrate in Sec. [TV[ 
If we let z = e^ M , the grand partition function of Eq. ( |I0D can be written as 

2 = ^^/ = ^^'^^ (14) 

N N {n} 

where the prime in the sum over {n} implies a sum over all n , ni, n 2 , ... = 0, 1, 2, . . . such 
that J2 P n p — N. If we define x = e _/3A , the last exponential in Eq. (14) can be written as 
x M with M — E/A — J2 P P n p and Z becomes 

z = E E = E E MM) x m z n , (is) 

N {n} N M 

where c^(M) is the degeneracy factor for the energy state M of N particles. 

In the ID oscillator problem, this degeneracy factor cn(M) has a very interesting mathe- 
matical property.El'§Si3ii It is just the number of ways that one can partition the integer 
M into N or less integers. For example, M = 4 can be partitioned in 5 ways: 1 + 1 + 1 + 1, 
1 + 1 + 2, 1 + 3, 2 + 2, and 4, which is equivalent to the number of ways that four Bose particles 
can be put in equally spaced states 0, 1, 2, 3, 4 to have 4 units of energy. Euler, Gauss, 
Hardy, Ramanujan, and many other famous mathematicians have contributed theorems on 
partitions.^ 

We do not need an explicit expression for cn(M). We can identify the canonical partition 
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function in Eq. (|15|) as 

Z n = J2cn(M)x m . (16) 

M 

We also have from Eq. (|10D that 

1 111 

z i z ) - n 1 _ zxP - 1 _ z i _ zx 1 _ zx 2 " ( i7 ) 

We next replace z in Z(z) by xz so that 

Z{xz) = 1 1 - • • • = (1 - z)Z(s), (18) 
1 — zx 1 — zar 

or 

£ £ MM) = (1 - z) Y, E c n(M) x m z n . (19) 

2V M N M 

We equate equal powers of z to obtain 

E c*(M) * M+Ar = £ c,v(M) x M - £ Civ-i(M) x M (20) 

M M M 



so that 



and 



x M Z N = Z N - Z N - U (21) 

(22) 
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This recurrence relation is trivial to solve explictly. We obtain 

N i N 

ZN=\[—- k = X{Z i m- (23) 

k=l 1 x fc=l 

Many derivations and uses of this result are found in the current research 
literature.00@il0'il'ii although the result itself has been around for many years.!'i!S0izl 
The above derivation is from Ref. Equation ( p3|) is applicable to the ID harmonic Bose 
gas. However, in a recent article, Schonhammeiil showed that an almost identical rela- 
tion holds for the ID ideal harmonic Fermi gas. His derivation can easily be adapted to 
work for bosons. The Fermi partition function differs only in having a factor e _/3£;o , where 
Eq = N(N — l)A/2 is the Fermi zero-point energy. This result means that the internal 
energies, given by — din Z^/d(3, differ only by E and that the two systems have identical 
heat capacities@ii given by 

dE f (/?fcA)V fcA m) 



As shown in Fig. 2, this quantity is linear in T at low temperatures (that is, for A <C ksT <C 
To) and approaches Nks at T ^> To. The relation between fermions and bosons for this 
system was first pointed out in Ref . . 

There is a curious aside to this relation between Bose and Fermi partition 
functions.iil@i3@ The ID harmonic gas has a constant density of states, which leads to 
the equality of the heat capacities. Consider instead free particles in 2D where the single- 
particle states are e p = p 2 /2m and the density of states is a constant because pdp = mde. 
One can show by using standard grand canonical ensemble techniques that the 2D free Fermi 
and Bose gases have identical heat capacities. This result is rather remarkable considering 
the considerable difference between the Fermi and Bose derivations of C. Moreover, the 2D 
Fermi/Bose heat capacity is fit extremely well by Eq. (f24|) . 



IV. MORE RECURRENCE RELATIONS 



It is possible to go beyond the ID harmonic case and derive canonical recurrence relations 
valid for any ideal gas. We first derive relations for the distribution functions following 
a method due to Schmidt, El who showed how the standard Fermi and Bose distribution 
functions in the grand canonical ensemble could be obtained by this means. We have 

n p (N) Z N = E^ e " /3Efc£fcnfc ^An I , (25) 
W 

where the Kronecker delta restricts the sum to N particles. Let n p = n' p + 1 and n' k = 
for k p. Then 

n p (N)Z N = + l)e-^-<+^_ liSX . (26) 

{«'} 

The term corresponding to n' p = — 1 does not contribute and the right side involves standard 
partition function sums corresponding to N — 1 particles. We have 

n p (N) Z N = [n p (N - l)Z N _ x + Z N _ X ] e"^, (27) 

or 

n p {N) = e -^%± [1 + n p (N - 1)] , (28) 

which is Schmidt's recurrence relation.^ Equation (p8|) appeared much earlier in the 
literature.00El@ These derivations assume that n p (N — 1) w n p (N) and use the rela- 
tion Z N = e l3FN , where F N is the Helmholtz free energy and F N — F N _i w dF N /dN = /i, to 
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find 

n ' {N) = ef*-,zjz^ - 1 = SJE^T < 29 > 

We end up with a canonical derivation of the standard grand canonical ensemble distribution 
function for bosons. An analogous derivation is valid for fermions. 

Unfortunately, the assumption that n p (N — 1) rs n p (N) leads back to the same fluctu- 
ation inaccuracies inherent in the grand canonical ensemble. But we need not make this 
assumption; Eq. fl28|) has a direct solution. By using the obvious starting values, n p (0) = 
and Zq = 1, we can prove by induction^! that 

N 7 

n p (N) = Y,e^ l ^p. (30) 
i=i Z N 

To use this relation, we need the partition functions involved. If we sum the relation over 
all p, we get 



v 



E 



(31) 



The quantity in square brackets is just the one-body canonical partition function at the 
effective inverse temperature (51, and a recurrence relation for results: 

1 N 

Zn = ^Y. Z ^ 1 ) Z n-i- (32) 
iV i=i 

This relation was apparently first derived by Landsberg,@ but appears many times in the 
current research literature.ilii0'0'0Ei@0'0'0iiii It was derived in this journal by Ford,@ 
although he made no application of it. Ford showed that such a result stems from the 
relation of Fermi and Bose partition functions to symmetric polynomials. Recently Schmidt 
and SchnackS extended this idea. 

The use of Eqs. (^8|) and (|32| ) allows us to determine the canonical distribution functions 
for the finite ID harmonic Bose system. (We could as easily find the properties of the finite 
3D harmonic Bose gas.) We start the recurrence in Eq. fl3"2| ) with Z = 1 and find every 
Zl,L < N. These values are then put into (P5|), starting with n p (0) = to find each 
distribution function n p (L), L < N, in sequence. The results are shown in Fig. 3. There is a 
small disagreement between the results for the canonical ensemble and the grand canonical 
ensemble for iV = 500, but these become smaller for larger N. The real difference between 
the two ensembles arises in the fluctuations. 
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We can, in the same way, develop a recurrence relation for the mean square distribution. 
We find 

rf(N) = e -**%=i [l + n p (N - 1) + n p (N - 1) 2 1 . (33) 

With this relation and Eq. (|28|) we obtain the root-mean-square fluctuation in the ground- 
state distribution function for the canonical ensemble, which can be compared to the same 
result for the grand canonical ensemble as given by Eq. (p|) . The results are shown in Fig. 4. 
As expected, we see that the grand canonical ensemble result goes to the total number of 
particles as T becomes small, while that of the canonical ensemble goes to zero after peaking 
near the quasi-transition temperature To. 

We can sharpen the distinction concerning fluctuations by deriving the probablility Po(n) 
of finding n particles in the state with eo = 0. For the ID harmonic gas in the canonical 
ensemble this quantity is given by 



1 



Z N 



no,ni, 



1 J2 e^ A ^ +2 ^+3^ + -) e -/3A(n 1+ n 2+ n 3+ ...) 57v _ n ^ inr (34) 

Zn ni ,n 2 ,... 



The series in the second exponential on the right can be replaced by N — n and extracted 
from the sum; the remaining factor becomes a partition function for N — n particles or 

P ( n ) =e -/3A(v-n)^v-n_ (Canonical ensemble) (35) 
Zn 

The equivalent quantity can be derived for the grand canonical ensemble. We find 

P ( n ) = (1 - e^)e nl3fl . (Grand canonical ensemble) (36) 

The two Po{n) functions are plotted in Fig. 5 for T = 0.12T . We see that the canonical en- 
semble value is sharply peaked around the average value while the grand canonical ensemble 
function is incorrectly wide and monotonic. 

If we include interactions, the fluctuations in the grand canonical ensemble are tem- 
pered to a more physically reasonable value.0 However, we would still prefer to be able 
to treat the ideal gas correctly for its conceptual and pedagogical importance and because 
experimentalists can now reduce the effective interactions to near zero by use of Feshbach 
resonances.! 
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One might argue that the micro-canonical ensemble, in which the system of interest is 
isolated, is more nearly equivalent to the actual experimental conditions of Ref. [I] than the 
canonical ensemble. However, it can be shown by use of the recurrence relations for the 
microcanonical ensemble©!^ that the microcanonical and canonical ensembles give almost 
identical results. 



V. PERMUTATION CYCLES 

The sums in Eqs. (|j) and (|32D have a physical interpretation as sums over permutation 
cycles. This view of BEC was first developed by Matsubarall and Feynman,! and was 
recently discussed in this journal by one of us.@ Here we examine permutation cycles in 
the context of recurrence relations. The boson A-body wave function is symmetric and 
the density matrix can be written in terms of symmetrical permutations of particles.i The 
partition function is the trace of the density matrix and involves a sum over all permutations, 

Zn = ivT ^ J dr i--- dr N (r P1 , . . . , r PJ v|e _/9H |ri, ...,r N ), (37) 

where the variable rp j represents the coordinate of the particle interchanged with particle j 
in permutation P. 

Any A-particle permutation can be broken up into smaller permutation cycles.iJ^ For 
example, for N = 7 we might have a 3-particle permutation cycle 1 — > 2 — > 3 — ^ 1 plus a 
4-particle cycle 4 — > 5 — > 6 — > 7 — > 4. The corresponding matrix element in Eq. ( |3"TD breaks 
up into a product of cycle terms: 

(3125674|e- /3// |1234567) = (312 1 e"^ 3 1 123) <5674| e"^ 4 14567> 

= £ e' 3 ^ ]T e" 4 ^ = Z 1 (3(3)Z l (4(3) } (38) 

p m 

where, for example, H 3 = hi + h 2 + ^3 with each hi being a one-body Hamiltonian. We have 
reduced the cycle matrix elements to one-body partition functions at an effective tempera- 
ture. The details of this derivation are given in Ref. |59|. Every term in the sum of Eq. ( |37j) 



can be reduced in this way to a product of permutation cycles represented by products of 
one-body partition functions. A single configuration will consist of q\ loops of length 1, q 2 
loops of length 2, etc., and may be arranged in C(qi,q 2 , ■ ■ ■) different ways. Thus we can 
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write 

1 

ivT 



n = — E C(gi,g 2 ,...)II^WS (39) 



{91,92,—} ' 

where the prime on the sum implies that it is over all combinations of permutation cycles 
such that 

5>Z = JV. (40) 
i 

Feynmani has given an argument (repeated in Ref. |59|) to show that 

°(«>».-)= l « v .. N l\»\..: (41) 

Again there is a connection to the theory of numbers. The sum in Eq. ( |39| ) is over the number 
of ways of partitioning the integer N into smaller integers. An example of C(qi, (fe, . . .) is 
the breaking of particles 1, 2, . . . , 5 into a 2-cycle and a 3-cycle, that is, partitioning 5 into 
2 + 3. With five particles there are several ways of doing this: We can take particles 1 and 
2 in the 2-cycle with 3, 4, and 5 in the 3-cycle, or take particles 1 and 3 in the 2-cycle with 
the remaining particles in the 3-cycle, etc. In all there are C = 5!/(2 1 3 1 l!l!) = 20 distinct 
ways of doing this, as the reader can verify. 

For the case of the ID harmonic Bose gas, we can combine Eqs. (f23"l), and ([39"1) to 
find an interesting relation: 

N I 

T\Z 1 (Bk)= y' —, TT^iW (42) 

k=l {91,92,—} H H 1 



J N — 

k=l 



or 



(1 - x)(l - x 2 ) • • • (1 - x») parti £ s ofN 1*2* ... gi \q 2 \...(l- x)*(l - x 2 )^ • • • ' (43) 

Equation ( ^3|) is known as Cayley's decomposition in the theory of partitions of numbers.il 
A simple special case is 1/(1 — x)(l — x 2 ) = | [1/(1 — x) 2 + 1/(1 — x 2 )]. 

The right-hand side of Eq. ( ff2]) is, in fact, the solution of the recurrence relation (p2|). 
This solution tells us how to interpret the recurrence relation itself. The sum in Eq. fl52"|) 
is a sum over permutation cycles: We can generate the iV-body partition function for the 
Bose system by adding a particle either as a singlet (Zi((3)) with the other N — 1 particles 
grouped independently (Zjv-i), or as part of a pair-exchange cycle (Zi(2/3)) with the other 
N — 2 particles in all their possible combinations (Zjv_2), or as part of a triple cycle, and so 
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on, with each configuration having equal probablity 1/N. (The work of Laloe et al. on 
interacting gases is closely related to this approach.) 

A further useful quantity is the average number p { = q x l of particles involved in permu- 
tation cycles of length I. This number is found by using Eq. ([32]): 

N = T,Pi = T,Zm)^ (44) 
i i Z N 

which tells us that 

V l = Z 1 (pi)?£=L. (45) 

We can easily plot p~ l , but before we do, it is useful to separate out the contributions to 
p x from the condensate and the excited states. From Eq. ( p0[) we have 



-o = E^ (46) 

l Z N 



so that pi = Zn-i/Zn is the contribution of the condensate to the average particle number 
in the permutation cycle of length /. However, this quantity is essentially unity until / is of 
order N. On the other hand, for small I in the one-dimensional oscillator, Zi(/3l) ~ l/(/3/A) 
in Eq. ( [I5| ) (cf. Eq. (|^)), which is the contribution of the non-condensate to the permutation 
cycles. Figure 6 shows p x vs. I for the canonical ensemble. Notice the rapid drop-off of p { for 
small I, corresponding to the non-condensate. However, for low temperatures (that is, when 
there is a large condensate) p x ~ pj ~ 1 out to a value equal to the condensate number n , 
where p x must drop off to satisfy Eq. fl46|) . We might have guessed before the calculation that 
the condensate consists only of very long permutation cycles of approximately n particles. 
Now we see that this is not true; the condensate particles have equal probability of being in 
singles, pair cycles, triple cycles, and so on out to an n -cycle. 

The dotted line in Fig. 6 is the grand canonical ensemble estimate of p t . We find this 
estimate from Eq. (Q), which we can showi! also to be a sum over permutation cycles. Thus 

Pi = e ll3fJl Zi(f3l). (Grand canonical ensemble) (47) 

We also have 

J oo 



n 



— - = e 13 ^ 1 . (Grand canonical ensemble) (48) 



e ■*• i=i 



The summand is the condensate contribution to p z in the grand canonical ensemble. 
The dotted line in the plot shows that the grand canonical ensemble does not do a very 
good job of representing the true nature of condensate permutation cycles in the Bose gas. 
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Because the condensate contribution to p L must drop off at n , we could estimate n by 
finding the value of I for which p l « 0.5. For the lowest temperature 0.12To in Fig. 6 this 
estimate gives no = 472, whereas the exact result is 471. This approach also allows us to 
estimate n in path-integral Monte Carlo simulations involving trapped interacting particles, 
where no standard estimators of no exist .3iE 

The picture we have then of the condensate is that it does indeed fluctuate wildly, not in 
overall particle number as would be suspected from the grand canonical ensemble, but rather 
in how it breaks into permutation cycles, with its particles having an equal probability of 
being in cycles of all sizes up to the condensate number itself. 



VI. DISCUSSION 



Our goals in this paper have been multiple: (a) Illustrate the inadequacies of the grand 
canonical ensemble in its depiction of fluctuations in a Bose-condensed system, (b) Show how 
the canonical ensemble can do a much better job in describing a condensed system, (c) Find 
recurrence relations that allow simple treatments of finite ideal Bose systems, (d) Study 
a simple model system, the ID harmonically trapped ideal Bose gas, which illustrates all 
of the important elements of BEC in a finite system and has its own particularly simple 
recurrence relation, (e) Delve deeper into the intricacies of BEC to find the physical meaning 
of the recurrence relations by looking at permutation cycles, (f) Illustrate the close relation 
between Bose and Fermi systems for cases where the density of states is constant. 

All our work here involves non-interacting Bose systems. Although real gases have non- 
negligible physical effects due to interactions, many of the ideas we have developed carry over 
to the interacting regime. 0E3E3EB0E3Q Although it is probably not possible to develop 
recurrence relations for interacting systems,El the idea of the condensate involving all sizes 
of permutation cycles, the usefulness of even a ID model of a Bose gas, and the close relation 
between the ID Bose and Fermi gases, are ideas that are still expected to hold for interacting 
systems. 

An important feature of the present paper is that the ID model is so simple that instruc- 
tors can use it in elementary courses in statistical physics without sacrificing much important 
physics. The computer programs needed to carry out the recurrence relations are simple 
and can be coded by the students themselves. There are thought to be few models where 
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the canonical ensemble is soluble, but we have seen here that any ideal system where the 
single-particle energies are known is actually tractable. Usually one does not care so much 
about using the canonical ensemble because the grand canonical ensemble makes the math 
easier. However, we have seen that for the case of BEC the grand canonical ensemble is not 
always accurate and the canonical ensemble becomes not only accessible but necessary. 
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Mullin Fig. 1 

FIG. 1: Grand canonical ensemble calculation of the number of particles in the two lowest states 
versus T/Tq for the ID harmonic Bose gas. The results shown for all the figures are for iV = 500. 
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Mullin Fig. 2 

FIG. 2: Heat capacity per particle (in units of fee) for the ID ideal harmonic Bose gas. This 
quantity is identical to the same quantity for a ID ideal harmonic Fermi gas.ll 
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Mullin Fig. 3 



FIG. 3: Comparison of the grand canonical ensemble and canonical ensemble calculations of the 
number of particles in the two lowest states versus T/Tq for the ID harmonic Bose gas. 
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Mullin Fig. 4 

FIG. 4: Root-mean-square fluctuation of the number of particles in the ground state of the ID 
harmonic Bose gas for the grand canonical and canonical ensembles. The fluctuations in the 
condensate in the grand canonical ensemble become as large as the occupation itself, which is 
unphysical. The canonical ensemble result is more reasonable. 
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Mullin Fig. 5 

FIG. 5: The probability of finding n particles in the ground state versus n for the ID harmonic 
Bose gas at T = 0.12To for both the canonical ensemble (solid line) and grand canonical ensemble 
(dotted line). The result for the grand canonical ensemble is unphysical. 
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Mullin Fig. 6 

FIG. 6: The number of particles in permutation cycles of length I versus I for the canonical ensemble 
(solid lines) at various temperatures for the ID harmonic Bose gas. Also shown by the dotted line 
is the same quantity at the lowest temperature for the grand canonical ensemble. 
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